Enhanced damping of ion acoustic waves in dense plasmas 
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A theory for the ion acoustic wave damping in dense plasmas and warm dense matter, accounting 
for the Umklapp process, is presented. A higher decay rate compared to the prediction from the 
Landau damping theory is predicted for high-Z dense plasmas where the electron density ranges 
from 10 21 to 10 24 cm and the electron temperature is moderately higher than the Fermi energy. 
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The ion acoustic wave, a longitudinal collective mode 
in plasmas, plays a crucial role in a range of applica- 
tions, such as the Thomson scattering [3, [J and the Bril- 
louin scattering Q. Understanding the dynamics of the 
wave in dense plasmas or warm dense matter is impor- 
tant in various context, including the inertial confinement 
fusion 0, [I[ and the compression of x-rays [fHU . 

The decay rate of the ion acoustic waves in plasmas is 
often modeled by the prevalent Landau damping theory. 
However, this theory is inadequate for dense plasmas, as 
the Umklapp process, which is not accounted for, be- 
comes pronounced in high densities. It was shown that 
the Umklapp process dominates the Landau damping for 
low-fc plasmons [Tol - [T2j . Even though the detailed under- 
lying physical mechanisms of the plasmons are different 
from those of the ion acoustic waves, it is expected that 
the Umklapp process is also important to the ion acoustic 
waves. The goal of this paper is to estimate the effect of 
this process. Starting from the plasmon damping theory 
in dense plasmas [11], a new theory predicting the ion 
acoustic wave sampling is proposed and a regime where 
the decay rate is larger than the prediction by the Lan- 
dau damping theory is identified. This result would have 
implications on the Brillouin scattering of dense plasmas, 
the x-ray Thomson scattering, and the reflection problem 
in the inertial confinement fusion. 

First we provide a brief review on the Landau damping 
theory for the ion-acoustic wave. Only a neutral plasma 
of single ion-species ions is considered for simplicity. We 
denote the electron (ion) temperature by T e (Tj), the 
corresponding density by n e (fij), the mass by m e (roj), 
and the charge by Z e = 1 (Zi = Z), where the charge 
neutrality condition reads n e = Zrii. The longitudinal 
dielectric function is e(k, uj) = 1 + \e + Xii where 
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frequency, and fc >e is normalized as J f^ e = lc? 3 v. We 

assume va < cj/k <C v te , where Hti,tc = \/T\tJwH^e- This 
is a necessary condition for a moderate Landau damp- 
ing. Under this assumption, \e and Xi can be esti- 
mated to be Xe — l/(^Ade) 2 and Xi — — {ui/u) 2 , where 
Ado = \/Tej '4:Tm e e 2 is the Debye screening length. The 
condition e = yields the dispersion relation for the ion 
acoustic wave, 
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where V = y/ ' ZT e /rrii. The ion acoustic wave decay rate 
from the Landau damping theory is given to be 
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and S = w iaw /fcw t o,ti- 

According to the Landau damping theory for the ion 
acoustic wave, there are many electrons satisfying the res- 
onance condition, of which energy is very low compared 
to the average kinetic energy. The distribution function 
around the resonance condition hardly varies, and nearly 
even electron population on both sides of the resonance 
condition results in little net energy transfer between the 
wave and the electrons. In other words, the derivative 
of the distribution function with respect to the velocity 
nearly vanishes at the resonance condition, which makes 
the decay rate small. This physical picture that electrons 
are freely-streaming and interacting only with the wave is 
no longer accurate in dense plasmas, because the distor- 
tion in the electron motion due to the presence of the ions 
(i.e., the Umklapp process) becomes important [ToL [TTI| . 

The effect of the Umklapp process on the ion acous- 
tic wave damping can be analyzed by an approach simi- 
lar to what was taken in developing the plasmon damp- 
ing theory for dense plasmas in Ref. [ljj]. Below we 
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follow nearly the same steps as in Sec. IV therein. 
In the presence of the potential of the form 0(x, t) — 
</>exp(ik-x — iujt) + <fi* exp(— ik-x + iwi), the wave packet 
is modified to be 
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where the perturbation theory of the first order is used 
assuming the perturbation is weak, | ) denotes the origi- 
nal eigenstate, |j ) denotes the perturbed eigenstate, and 
a is an index for the eigenstate. Then the perturbation 
in the density is given to be 

5n(k,u) =J2f(°i)\(°i II exp(ik-x) || a,)] 2 

= 2J e ^(/( CTl ) ~ /( Cr 2))/3(cri,(T2,k,w), 

where /(<7j) is the occupation number, and /3(cti,<72) is 
(<ti| exp(— zk • x)|<72)(cr 2 | exp(zk • x)|oi) 
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With an appropriate choice of the eigenstates under a 
given condition, Eq. (Q} can b e ap plied to various sit- 
uations. For example, Sturm used the eigenstate 

W) = |q) + E qi#q lqi)(qil^|q)/(^ q - £ qi ), where q 

and qi are the wave vectors, and (qi|V|q) is the pseudo- 
potential. We will suggest an appropriate eigenstate for 
the ion acoustic wave, after discussing the ion dynamics 
below. 

Assuming each ion is located at Xj in the time scale 
of the electron damping on the wave (i.e., the Born- 
Oppenheimer approximation), we compute the damping 
rate of the ion-acoustic wave due to the electrons in the 
presence of the spatially fixed ions, and then obtain the 
average dynamics by integrating over the probability dis- 
tribution of Xj. Without loss of generality, Xj can be 
assumed to follow the correlation average, 
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FIG. 1: The decay rates computed by the theory developed 
in this paper, compared with the prediction of the Landau 
damping theory (dot-dashed line) which is almost constant 
for a range of the wave vector. The parameters used are 
nii = 26, Z = 3, Ti = 20 eV, and T e = 100 eV. Two cases 
of the ion acoustic wave vectors are shown, k — 0.1 kdc and 
k — 0.4 fcdo- n e is in the unit of cm -3 . 



the ions. In the presence of spatially fixed ions, the elec- 
tron's free wave eigenfunction is modified to be 



k) = |q) + E 



rf 3 qi exp(-iqi • Xi)U (q - qi) 
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where U(q) = 47r Z e 2 / (|q| 2 + k 2 c ) is the Fourier transform 
of the screened ion-electron potential. 

For the plasmon damping, E e f / (q, ch) = E(q) — E(c[ 1 ) 
can be used, as was used by Sturm [l(|. This is a good 
approximation for the plasmons, as the plasmon energy 
huj e is very high and the non-degenerate perturbation 
theory is good enough for the calculation of the distor- 
tion of the electron wave packet due to the presence of 
ions. However, in the case of the ion acoustic waves, 
the wave energy ?tWi aw is much smaller than the electron 
temperature, and the perturbation theory is almost de- 
generate. In order to account for this effect, we choose 
#e//(q,qi) = E(q) + \E(q)-E(q 1 )\, which is a good ap- 
proximation for the nearly elastic or large inelastic scat- 
tering. The susceptibility obtained from Eqs. (HJ, ([5]) and 
© is 

Xe(k,w) = Xrpa(k,w) +Xdonsc(k,w), 

where Xrpa is given in Eq. (JTJ) , and the subscript stands 
for the random phase approximation. Xdense is of our 
main interest, which is given to be 

Xdc„se(k, W ) = —g- J pnyM*) 

d 3 d U 2 (s) /(q + k + s)-/(q) 

(2tt) 3 A 2 (q, k, s) hw - E{q + k + s) + E(q) ' 
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(5) where A(q, k, s) is 



where V c is the volume of the region under consideration, ^ (l) kj s ) 
and n/(s) is the static two-point correlation function of 
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The imaginary part of Xdcnsc can be obtained from 
Eq. by replacing the denominator, l/(frw~ E(c[ + k + 
s) + E(q)), by the delta function TrS(huj — E(q + k + s) + 
E(q)). For a Maxwellian plasma, Immense] can be fur- 
ther simplified using the velocity integration v = /iq/m e , 
instead of the wave vector q: 

f n z (s)d 3 s ( 4^e 2 V ^ 

Ini[XdcnscH i T^l^TfciJ If J 

x y rf 3 v (A" 2 (v, k, s)tt(5(w - (k + s) • v)) 

/m(v + ^)-/m(v-^) 
x if > l°J 

me 

where /a/ is the Maxwellian distribution of the electron 
temperature T e , satisfying the normalization condition 
/ /m^ 3 v = 1. A~ x (v, k, s) is given to be 

A 1(v ' k ' s) = S(v _k±s ) _| i?(v _k±s ) _ £ ; (v _k_s ) | 

1 

" E(v + ^) - |£(v + i=i) - £(v + ^)| 

where s can be integrated in the spherical coordinate sys- 
tem. For a given set of k and s, the integration over the 
velocity can be reduced to a two-dimensional integral, as 
one variable is eliminated by the delta function. The de- 
cay rates given by Eqs. (|3]) and (JSj) (Fig. [IJ exhibit that 
the newly computed decay rate is much higher than what 
is given by the Landau damping theory when fcAdo = 0.4. 
Similar integrations for various physical parameters show 
that the regime where the Umklapp process is important 
has the electron density ranging from 10 21 to 10 24 cm~ 3 . 



The high ion charge would even further enhance the de- 
cay rate. In the presence of high-Z ions, the electrons 
interact more strongly with the ion acoustic wave, the 
ions being used as the momentum storage. 

In summary, it is suggested that the ion acoustic wave 
decay rate could be much higher than the prediction by 
the Landau damping theory (see Eq. (JSJ). A rather 
rough theory, accounting for the effect of the Umklapp 
process on the decay, is presented. The theory proposed 
here is far from being complete. For instance, the per- 
turbation expansion given in Eq. © may deviate signifi- 
cantly for the electrons of small kinetic energy (less than 
1 eV). However, it highlights an important weak point of 
the prevalent Landau damping theory, when applied to 
dense plasmas: The ion acoustic wave damping in dense 
plasmas may arise from the three-party interaction (elec- 
tron, ion, and wave), as opposed to the two-party inter- 
action (electron and wave) the Landau damping theory 
is based on. 

Some comments on the relationship of our theory to 
the plasma kinetic theory are in order. In plasma physics, 
the damping of various waves due to the ion-electron col- 
lisions is often considered in terms of the drag of the 
electron motion due to ions. For example, the inverse 
bremsstrahlung is treated in this way hj. When the 
classical ion-electron collision picture breaks down, the 
dielectric function approach presented here is superior 
to the classical kinetic approach as demonstrated in the 
computation of the plasmon damping in the solid state 
physics problems [Hf, El- In the regime we consider, the 
classical electron-ion collisions are no longer valid and the 
quantum mechanical treatment is necessary. It would be 
interesting to consider the quantum ion acoustic wave as 
well [HI], in the context of the quantum electron degen- 
eracy and diffraction [T34l8| . 



[1] S. H. Glenzer, L. M. Divol, R. L. Berger, C. Geddes, 
K. R. Kirkwood, J. D. Moody, E. A. Williams, and 
P. E Young, Phys. Rev. Lett. 86, 2565 (2001). 

[2] J. Chihara, J. Phys: Condens. Matter 12, 231 (2000). 

[3] J. R Murray, J. Goldhar, D. Eimerl, and A. Szoke, IEEE 
J. Quantum Eletronics 15, 342 (1979). 

[4] M. Tabak, J. Hammer, M. Glinsky, W. Kruer, S. Wilks, 
J. Woodworth, E. Campbell, M. Perry, and R. Mason, 
Physics of Plasmas 1, 1626 (1994). 

[5] R. P Drake, and S. H. Batha, Phys. Fluids B 3, 2936 
(1991). 

[6] M. Hentschel, R. Kienberger, Ch. Spielmann, G. A. Rei- 
der, N. Milosevic, T. Brabec, S. Corkum, U. Heinzmann, 
M. Drescher, and F. Krausz, Nature 414, 509 (2001). 

[7] P. Emma, K. Bane, M. Cornacchia, Z. Huang, H. Schlarb, 
G. Stupakov, and D. Walz, Phys. Rev. Lett. 92, 074801 
(2004). 

[8] V. M. Malkin and N. J. Fisch, Phys. Rev. Lett. 99, 
205001 (2007). 



[9] V. M. Malkin, N. J. Fisch, and J. S. Wurtele, 

Phys. Rev. E. 75, 026404 (2007). 
[10] K. Sturm, Z. Physik B 25, 247 (1976). 
[11] S. Son, S. Ku, and S. J. Moon, Theory of plasmon decay 

in dense plasmas and warm dense matter Phys. Plamas 

In press, (2010). 
[12] W. Ku, and A. G. Eguiluz, Phys. Rev. Lett. 82, 2350 

(1999). 

[13] J. M. Pitarke, and I. Campillo, Nucl. Inst. Meth. Phys. B 

164-165, 147 (2000). 
[14] J. Dawson, and C. Oberman, Phys. Fluids 5, 517 (1962) 
[15] F. Haas, L. G. Garcia, J. Goedert, and G. Manfredi, 

Phys. Plasmas 10, 3858 (2003). 
[16] S. Son, and N. J. Fisch, Physics Reviews Letters 95, 

225002 (2005). 

[17] S. Son, and N. Fisch, Physics Letters A 329, 76 (2004). 
[18] S. Son, and S. Ku, Physics of Plasmas 17, 010703 (2010). 



